Multiple-copy distillation and purification of phase diffused squeezed states 
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We provide a detailed theoretical analysis of multiple copy purification and distillation protocols 
for phase diffused squeezed states of light. The standard iterative distillation protocol is generalized 
to a collective purification of an arbitrary number of N copies. We also derive a semi-analytical 
expression for the asymptotic limit of the iterative distillation and purification protocol and discuss 
its properties. 



I. INTRODUCTION 



Continuous variable (CV) quantum information pro- 
cessing [l[ is complementary to discrete variable quantum 
information based on qubits. It employs the encoding of 
information into field modes of light 0, d, 0] or into the 
collective spin state of many particles such as in a cloud 
of atoms [1, H, 0| ■ Encoding quantum information into 
continuous variables of optical modes exhibits several ad- 
vantages over approaches based on single-photon qubits. 
For example, it allows for the deterministic realization of 
several important protocols, such as the generation of en- 
tanglement and entanglement swap ping [l2j], 
quantum teleportation I2J, |3|, quantum cloning [13J, and 
quantum dense coding [14j |. using only linear optics, op- 
tical parametric oscillators as squeezed light sources and 
balanced homodyne detection. It is also possible to store 
continuous variables of light in atomic memories, either 
by means of direct interaction [|| or by teleportation Q . 

In the description of CV systems, an important role 
is played by probability quasidistributions, which repre- 
sent each field mode by a real two-dimensional function 
[li|. When this function is of Gaussian shape, the cor- 
responding state is referred to as Gaussian. Due to the 
ease of both experimental and theoretical treatment of 
such states, they are, together with Gaussian operations 
preserving the Gaussian nature of the quasidistributions, 
a significant part of CV quantum information. How- 
ever, attempts to estabilish entanglement between two 
distant parties, which is one of the fundamental tasks 
quantum information is facing, are hampered by deco- 
herence and thus, purification and distillation of entan- 
glement are needed to remove the added noise and re- 
store, at least approximately, the initial entangled state 
[HI, [13] • Indeed, the entanglement distillation and pu- 
rification techniques allow to extract from many copies 
of weakly entangled mixed state a single copy of highly 
entangled almost pure state by means of local quantum 
operations on each party's side and classical communica- 
tion between the two parties sharing the states. In gen- 
eral, distillation refers to process of increasing entangle- 
ment (or squeezing), while purification is associated with 
reduction of the state's mixedness. Within the frame- 
work of this paper both effects involve each other and 



the corresponding terms will be used interchangeably. 
Unfortunately, purification and distillation of entangle- 
ment are impossible to accomplish for Gaussian states 
by using feasible Gaussian operations alone [HI, [l|| H(| • 
Note, that a similar no-g o theorem also holds for single 



mode squeezed states [21 



Namely, it is impossible to use 
passive Gaussian operations, balanced homodyne detec- 
tion and feedforward to transform an arbitrary number of 
Gaussian states into a single state with squeezing better 
than the initial one, where the squeezing is characterized 
by the lowest eigenvalue of the covariance matrix. On the 
other hand, Gaussian operations are sufficient to improve 
purity of a squeezed state if loss of squeezing is accepted 



To allow for increasing the entanglement of Gaussian 
states by local means, one has to at one poi nt step out of 
the Gaussian domain, as was utilized in |23l . l24l. l25ll . The 
complete distillation protocol, proposed by 25l 2a . Wf\ , 
consists of two steps: de-Gaussification by a single pho- 
ton subtraction [2J, [28|, [2^, [3(1 [H| followed by Gaussi- 



fication by means of interference on beamsplitters, mea- 
surement and conditioning. However, if the initial states 
are already non-Gaussian, for example due to the effect 
of non-Gaussian noise, it is possible to distill and pu- 
rify them by employing only linear optics, homodyne 
detection and post-selection [H, [H, [3J, [35[. Recently, 
we have proposed [lH and experimentally demonstrated 
[33l . HH distillation of squeezing from states disturbed 
by phase diffusion noise. This noise, caused by random 
fluctuations of optical phase, commonly occurs in optical 
communication links. The distillation protocol is based 
on the interference of two copies of the state at a bal- 
anced beamsplitter, followed by homodyne detection of 
one mode triggering the acceptance or rejection of the 
other mode that contains the purified state [25|, [26|, [52] • 
The procedure can be straightforwardly extended to ac- 
commodate two-mode squeezed states and can also serve 
as the Gaussification part of the universal CV distilla- 
tion protocol based on dc-Gaussification and subsequent 
Gaussification. 

In this paper we provide detailed theoretical numeri- 
cal treatment of the distillation and purification of phase 
diffused single-mode squeezed states. We consider the 
standard iterative purification protocol [32| as well as 
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FIG. 1: (Color online) Setup for purification of single-mode 
squeezed states. Two single-mode squeezed states are sent 
through independent dephasing channels which result in their 
decoherence. They are then mixed on a balanced beamsplit- 
ter (BS) and homodyned in both output ports. The virtual 
beamsplitters (VBS) are meant to simulate inefficient homo- 
dyning. 



the simultaneous collective purification of an arbitrary 
number of copies of the decohered state. We also dis- 
cuss the possibility of conditioning on the measurement 
of an arbitrary quadrature operator, up to abandoning 
the concept of definite quadrature completely and using 
phase-randomized homodyne detection. Finally we de- 
rive a semi-analytical expression for the state obtained 
in the limit of infinite number of iterations of the distil- 
lation procedure and discuss its properties. 



II. ITERATIVE PURIFICATION 

The single step of the purification protocol for single- 
mode squeezed states is depicted in Fig. [T] It consists 
of interference on a balanced beamsplitter (BS) followed 
by homodyne measurement of a single mode, detecting a 
value of the quadrature operator that was, prior to effects 
of the noise, squeezed [32], H^l- If the absolute value of 
the measured quadrature x is lower than a certain pre- 
selected threshold X then the purification succeeded and 
otherwise it failed. The purified state is present in the 
other output port of the beamsplitter and is available 
for further applications or for the next iteration of the 
purification procedure, whereby two purified copies of the 
state are purified again. The success of the procedure can 
be verified by suitable homodyne measurements on the 
purified state, as schematically shown in Fig. 1. We will 
treat the purification in the Fock basis which allows us 
to calculate the full density matrix of the purified state 
after k iterations of the purification protocol, 
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squeezed vacuum state can be expressed as 

P (a) = [ U(<t>)p G U\4>)$((t>)d(t>. (2) 



Here p G denotes density matrix of the initial Gaus- 
sian state, denotes the probability distribution of 
the random phase shift <fi and U (4>) is a unitary phase- 
shift operator with matrix elements (m\U\n) = e m ^5 mn , 
where the r5 m „ stands for the Kronecker delta. We can 
rewrite Eq. @ as follows, 
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where pG,m,n = (m\p G \n) and 
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(3) 



(4) 



is a Fourier transformation of $(</>). Throughout this pa- 
per we assume that the phase fluctuations have Gaussian 
distribution, 
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(5) 



On inserting this probability distribution into Eq. (TJJ wc 
obtain /„ = e^^l" 1 . 

The Wigner function corresponding to the squeezed 
vacuum state pa reads 



W(x,p) 



1 



(6) 



where V x and V p denote the variances of the x and p 
quadratures, respectively. For a vacuum state one has 



V x = V p 



and the state is squeezed when V x < -~ or 



V p < Tj. Note that due to the proper normalisation, all 



variables used throughout the paper are dimensionless. 
Presently, these values could go as far as V x « 0.06 - 
0.05 (~QdB to -10dB), V p m 16 - 25 (15dB - 17dB) 
[lElHl, but values around V x m 0.16 - 0.22 (-5dB to 
—3dB) are common in contemporary experiments (33l . 
HU, [35|. The density matrix elements PG,m,n can be 
evaluated by noting that the Husimi Q-function defined 
as Q(a) = -(a|/o|a), where \a) denotes a coherent state 
with amplitude a, is a generating function of the density 
matrix elements in the Fock state basis, 
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(7) 



The Q-function corresponding to the Wigner function © 
is also Gaussian, 



where |n) denotes the n-photon Fock state of the field Q( a > 
mode. The initial density matrix of a phase-diffused 



1 



■.exp[(U-l)aa*-T(c 



% (8) 



7T\ V x V p 



3 






0.8 








0.6 


■ — 




W 




W 


0.4 


CD 


O 




on 


0.2 


w 










0.2 


0.4 


0.6 

a 


0.8 


1 










(c) . 














0.2 


0.4 


0.6 


0.8 


1 



FIG. 2: (Color online) Performance of the four-step iterative 
purification protocol. The variance (a) and purity (b) of the 
purified state as well as the total success rate of the purifi- 
cation (c) are plotted as functions of the phase noise a. The 
curve for the initial phase diffused state (blue dashed line) is 
plotted along curves of several iterations (green solid lines), 
where the arrows indicate the direction of increasing number 
of steps. The parameters were r\ = 0.85, 14 = 0.2, V v — 2, 
and X = 0.45. 



where V x = V x + 1/2, V p = V p + 1/2 and 
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On inserting function ([SJ into expression ([7]) we obtain 
after some algebra 
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if m + n is even and PG,m,n = if m + n is odd. 

The distillation and purification scheme shown in 
Fig. 1 produces with certain probability from two copies 
of the state p^ -1 ' a single copy of a purified state p( k > 
and we can formally write the purification map as 



,(A0 



(11) 



As the formula suggests, the protocol can be iterated and 
the outputs of each step can be used as inputs for the next 



iteration. As we shall show, each iteration increases the 
squeezing and purity of the state and Gaussifies it. 

The interference of two copies of the state p^ 1 ) on 
a balanced beamsplitter is governed by a unitary trans- 
formation. Since this transformation preserves the total 
photon number, we can write it in the Fock basis as fol- 
lows, 



\mi,m 2 ) -> ]T]^m 1 ,m 2 l m i +a,m 2 - a), 



(12) 



where 
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dl(d + a)\(mi - d)\(m 2 - d- a)V 



(13) 



After mixing on a BS the first output mode is measured 
in a balanced homodyne detector where it is projected 
on the eigenstate of the x quadrature \x). The positive 
operator valued measure (POVM) element corresponding 

to conditioning on |x| < X reads C = J_ x \x)(x\ and its 
matrix elements in the Fock basis can be expressed as 



(m\x)(x\n)dx. 



(14) 



x 



This can be easily evaluated numerically by recalling that 
, , , 1 



y/wi2 n n\ 



H n (x)e~ x l\ 



(15) 



where H n (x) denotes the Hcrmite polynomial. For any 
realistic detector with efficiency 77 < 1 we can use the 
model where the inefficient detector is replaced by a 
beamsplitter with transmittance rj whose auxiliary input 
port is in the vacuum state and which is followed by an 
ideal perfect detector [lj|. This yields 



C m , n (v) = ^ B m<a (r))B nia (r))C m - a . 



(16) 



where 




B m , a = il ("; )77 (m - a)/2 (l-r7)°/ 2 . (17) 



We are now in a position to combine all the above ex- 
pressions and write down the purification map (jll| in 
the Fock state basis, 

(fe) ^ V V V n { - k ~ V > rfi-V A a A b 
r / j / j / y Fmi ,m rma,n 2 mi .ma-^ni ,wa 

mi, ni m 2 ,Ti2 a,b 

xC mi+a , ni +b{ri) \ m 2 - a)(n 2 - b\. 



(18) 
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FIG. 3: (Color online) Trade-off between the success rate of 
the four-step purification protocol and the resulting squeezing 
(a) and purity (b) of the state. The parameters were r] = 0.85, 
V x = 0.2,Vp = 2, and a = 0.5. 
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FIG. 4: (Color online) Fidelity of the purified state with re- 
spect to a Gaussian state with the same covariance matrix 
(green solid lines), the arrow indicates increasing number of 
iteration steps. The lowermost curve (blue dashed line) cor- 
responds to the initial de-phased state. The parameters were 
r) = 0.85, V x = 0.2, V p = 2, and X = 0.45. 



The density matrix (JT8J) is not normalized and its trace 
is equal to the probability of success (success rate) of k 
iterations of the purification protocol, 



P(k) = Tr\pW]. 



(19) 



Typical numerical results are shown in Fig. We can 
see that each step of the iterative purification decreases 
the variance of x, i.e. increases the squeezing. Also, 



the purity of the state V = Tr[p 2 ]/(Tr[p]) 2 increases af- 
ter each iteration of the protocol and approaches some 
asymptotic value, which is generally less than unity, c.f. 
FiglDJb). So our protocol purifies the state but it does not 
generally distill perfectly pure states from the initially 
mixed states. One important parameter determining the 
practical feasibility and usefulness of the scheme is the 
probability of success of the protocol, which is plotted 
in Fig. [2fc). We can see that for the chosen acceptance 
window X = 0.45 the success rate P(N) is very high, 
the probability of success of each iteration is about 50% 
and the total probability of success for a protocol in- 
volving four iterations is still about 10%. The trade-off 
between the performance of the protocol and the success 
probability is illustrated in Fig. [3] which shows the fi- 
nal variance and purity of the state after four iterations 
as a function of the total success rate P. Higher pu- 
rity and stronger squeezing are achievable at the expense 
of reduced P. Our numerical calculations suggest that 
it is suitable to choose the threshold X ss \/Vx, which 
achieves good purification and squeezing enhancement at 
a reasonably high probability of success. 

The purification procedure also Gaussifies the state. 
To quantify this, we evaluate the fidelity of the state 
p( fc ) after k iterations of the protocol with the Gaussian 
state pa that has the same mean quadrature values and 
covariance matrix as the state p( k \ The fidelity of two 
mixed states is defined as follows, 



F = 



(20) 



and it holds that F = 1 if and only if pc = p^ ■ A value 
F < 1 is a clear signature of a non-Gaussian character 
of the state p^ k \ The results are shown in Fig. [5] which 
confirms that the present protocol indeed Gaussifies the 
state. After four iterations, we have F > 0.999 so the 
state is almost perfectly Gaussified. 

So far we have considered only conditioning on the 
measurement of the initially squeezed x quadrature. 
However, the method works even for a measurement of an 
arbitrary quadrature q(6) = xcos8 + psin9 [HI], where 
the conditioning is again of the form \q(6)\ < Q. This can 
be straightforwardly implemented into our Fock-state ba- 
sis formalism by re-defining the matrix elements of the 
POVM C(rj) as follows, 



C m , n {'n,6) = C m ,„(r?)e 



i(m — n)9 



(21) 



Note, however, that there are two prominent quadra- 
tures, x = q(0) and p = q(ir/2). Conditioning on one 
of these quadratures yields, for separate ranges of <r, the 
best improvement of the purified state. Furthermore, 
since improvement in the purified state can be seen no 
matter which quadrature q{9) is measured (HJ, it is possi- 
ble to consider also a phase-randomized homodyning. In 
this case, Eq. (|2"Tj) has to be averaged over the random 
phase 6 and only diagonal elements survive, 



(22) 




0.2 0.4 0.6 0.8 1 




0.2 0.4 0.6 0.8 1 

a 



FIG. 5: (Color online) Iterative purification using condition- 
ing on \p\ < Q (a) or phase-randomized conditioning (b), with 
the parameter settings 77 = 0.85, V x — 0.2, V p = 2, Q = 0.45. 
The arrows indicate increasing number of iterations (green 
solid lines). The uppermost curve (blue dashed line) corre- 
sponds to the initial phase diffused state. 

The results of numerical simulations of the iterative pu- 
rification are given in Fig. [5l Note that the advantage of 
conditioning on p grows with the number of iterations. 
For example, if a = 0.5 then after the first iteration the 
variance is more reduced by conditioning on |x| < Q than 
on \p\ < Q. However, for a four-step iterative procedure 
this is reversed and it is slightly better to condition on 
\p\ < Q rather than on |x| < Q. The phase-randomised 
purification can be seen as a compromise between condi- 
tioning on |x| < Q and \p\ < Q - for large (small) values 
of a it yields results better (worse) than conditioning on p 
(x) but worse (better) than conditioning on x (p). In fact, 
values obtained by phase-randomised purification are al- 
most identical to those produced by random switching 
between x and p detections. 



III. COLLECTIVE PURIFICATION 

The iterative purification scheme analyzed in the pre- 
vious section requires at least 2 k copies for k iterations 
of the protocol. Here we analyze a more general collec- 
tive purification scheme that works for any number N of 
available copies of the state. In particular, this protocol 
can be employed when three copies of the state are avail- 
able and it reduces to the previous scheme when N = 2 k . 
The proposed setup is shown in Fig. [6l The N copies 
of the phase-diffused squeezed state are combined on an 
array of N — 1 beamsplitters with transmittances tj and 
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FIG. 6: (Color online) Collective purification of N copies of 
the state. The beams are combined on an array of TV — 1 
beamsplitters BSj and all output modes except for the last 
one are monitored with balanced homodyne detectors BHDfe. 

reflectances rj, 1 2 + ?' 2 = 1. All output modes except for 
the last one are monitored by balanced homodyne detec- 
tors which measure the x quadrature of each mode. The 
purification is successful if the absolute value of each mea- 
surement outcome is below the corresponding threshold, 
|xj,out| < Xj. The choice of the thresholds would de- 
pend on the transmittances of the beam splitters. Here 
we shall consider ideal homodyne detectors with unit ef- 
ficiency, 77 = 1, and conditioning on |xj ut| = f°r which 
the protocol achieves the best performance. In practice, 
this would correspond to choosing sufficiently narrow ac- 
ceptance windows, i.e. very small thresholds Xj. 

Let us for a moment fix the random phase shifts <j)j of 
each mode. Then the variance of the quadrature Xj of 
the j-th mode reads 

Vj = V x cos 2 <j>j + V p sin 2 , (23) 

and in what follows we shall not explicitly display the 
dependence of Vj on <pj for the sake of brevity. The joint 
probability distribution of the x quadratures of all N in- 
put modes factors into the product of marginal Gaussian 
probability distributions, 

N x 2 

Pi^{x) = I[-=e-*%. (24) 

j=l \ /Z7TV 3 

In the Heisenbcrg picture, the output x quadratures af- 
ter interference on the array of beamsplitters can be ex- 
pressed as linear combinations of the input x quadra- 
tures, 

JV 

gj.out = y~]Uj,ixi, (25) 
1=1 

where Uj.i denote elements of a real orthogonal matrix 
which describes the action of the interferometer com- 
posed of the beamsplitters. For the scheme shown in 
Fig. Hit holds that 

JV-l 

U Ntl = n_r JJ ti. (26) 

3=1 
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Here we formally define tq = 1. The input-output rela- 
tion (|25[) can be inverted and we have 



N 



N 



(27) 



z=i 



The joint distribution of the output quadratures can be 
obtained by inserting the Xj expressed in terms of xi.out 
in Eq. |24|) . Taking into account that the Jacobian of 
the transformation (1231) is dct U = 1, we obtain 



JY 



-^joint (^out) J 



i 



cxp 



1 N 



j *£/,out ) 



1=1 



(28) 

The un- normalized distribution of Xjv,out conditional on 
projections xi, out = 0, I = 1, . . . , N — 1 then reads 



-F'cond( a; A r ,out) OC 

where 



1 



(29) 



V f-f 



N V 2 



i=i 



(30) 



The variance of a;jv i0 ut can be calculated by averaging 
properly weighted V over the N random phase shifts, 

1 r T/ 3 / 2 N 

v '" = ^iuL^U m) ^ <31) 



Here d(j> = d<j)i ■ 
stant, 



and A/jv is a normalization con- 
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v 1 ' 2 
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Ollli^ 



A/2 



(32) 



Two choices of the beamsplitting ratios are of particu- 
lar interest. Using balanced beamsplitters is experimen- 
tally most straightforward, since these beamsplitters are 
most commonly used in the lab. In this case we have 
tj = rj = 2- 1 ' 2 which implies U N>X = 2-^ N - 1 ^ 2 and 
U N ,j = 2~ {N ~ l+1 V 2 , I > 2. Notice that the output 
quadrature xjv.out is an unbalanced combination of the 
input quadratures. The second choice is to use unbal- 
anced beamsplitters which yield a balanced superposition 
of the input quadratures, Un,i = This happens 

if tj = + 1). This latter scheme becomes fully 

equivalent to the /c-step iterative purification procedure 
discussed in the previous section when N = 2 k . Note 
that in the iterative scheme we can also first let all the 
beams interfere on the balanced beamsplitters and then 
measure on all output ports except one. If we condition 
on xi ou t = then the two schemes become immediately 
fully equivalent. For realistic finite thresholds Xj the 




FIG. 7: (Color online) Variance of the x quadrature of the 
state obtained by purification from 2, 3, and 4 copies is plot- 
ted in dependence on a. The curve labeled "in" shows the 
variance for the initial state. The initial parameters were 
V x = 0.2, V p = 2, rj = 1 and X = 0. 



schemes become equivalent if conditioning on appropri- 
ate linear combinations of a^out is adopted. 

The results of numerical simulations of the collective 
purification procedure are shown in Fig. [7] We can see 
that each additional copy helps to further suppress the 
noise and reduce the quadrature fluctuations. The dif- 
ference between the performances of the two-copy and 
three-copy protocols is sufficiently large so that it should 
be possible to observe it in present-day experiments. 
The numerics also reveal that the difference between the 
schemes using balanced and unbalanced beamsplitters is 
very small and the scheme is quite robust and insensi- 
tive to changes in reflectances or transmittances of the 
beamsplitters of the order of a few percent, which is an 
important practical advantage. 

A conditional measurement in the form of a de- 
tection of arbitrary quadrature q(9) can be also con- 
sidered for the collective purification. To treat in 
a unified way conditioning on measurements of arbi- 
trary quadratures of the first N — 1 output modes, 
we shall present a more detailed treatment of the sce- 
nario, based on the formalism of covariance matrices 
(CM). Let £ = (xi,pi, X2,P2, ■ ■ ■ , Xn,Pn) denote the 
vector of quadrature operators. The covariance matrix 
Sj/c = ^({A£j, A£fc}) comprises the second moments. 
The CM of the single-mode squeezed state is diagonal, 
^SMS = diag(V^, V p ). A random phase shift <fi transforms 
the CM to S SM s(0) = i?(</>)S SM s-R T (</>), where 



R{4>) 



COS( 

sin t 



— sm ( 
cos 6 



(33) 



For fixed random phase shifts <j>j , the covariance matrix 
of the N-mode Gaussian state impinging on the array of 
N — I beamsplitters is given by 



JY 

0S SMS ( 



(34) 
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FIG. 8: (Color online) Collective purification from three 
copies of the state. The figure shows the variance of the x 
quadrature before (blue dashed line) and after purification. 
The quadratures measured by balanced homodyne detectors 
BHDi and BHD2 were the following: xi, 22 (green solid line), 
pi, pi (red dot-dashed line), xi, P2 (black dotted line). 

The covariance matrix of the output modes can be ob- 
tained from Ei n by a symplectic transformation. 

Sout = SBsSinSeS' (35) 

where the symplectic matrix Sbs governs the linear mix- 
ing of the quadrature components in the passive linear 
interferometer consisting of the array of N — 1 beam- 
splitters, sec Fig. |U To make our treatment simple, we 
shall again consider the limiting case of conditioning on 
cjj = 0, j = 1, . . . , N — L The quadratures which are 
being observed are specified by N — 1 relative phases 9j 
and we assume that all detectors have the same efficiency 
r\. We introduce a new covariance matrix 

Saux = vS(0)X out S T (6) + ^-=^J, (36) 

where / stands for the identity matrix, S(9) = 

©jLi an d we se t = 0. The single-mode phase 

shifts in the above formula map all measured quadra- 
tures onto the x quadratures. We want to calculate the 
variance of the quadrature xn conditional on all mea- 
surement outcomes being equal to 0. We first construct 
a submatrix of E a ux which comprises variances and co- 
variances of the x quadratures only, 

= S x 'S SiUX S x . (37) 

Here S x is a N X 2N matrix defined as Sj^j-i = 1 and all 
other elements arc equal to zero. The joint probability 
distribution of the N quadratures x = (xi, x 2 , . . . , xn) 
can be expressed as 

P(sb) = 1 ■ g-^s- 1 * (38) 

(2tt)^/ 2 ^E^ 

where the |.| denotes determinant of the matrix. The 
probability distribution of xn conditional on x± — 



X2 ■ ■■ = Xn-i = reads 

p ^ MK ^^ e ~*' (39) 

where Vn = ^nn] is the conditional variance of 

xn ■ Since we assume that the mean values of all quadra- 
tures are initially zero and the purification scheme pre- 
serves this property, the resulting variance of a; at after 
purification can be evaluated as properly weighted aver- 
age of Vn, 

Here d(f) = dtpi ■ ■ ■ d<f>N and A/jv is a normalization con- 
stant, 

*~=lMfk m)d *- <4i) 

Note that V ou t determined in this way corresponds to the 
variance of the purified beam that would be measured by 
balanced homodyne detector with efficiency r\. The ex- 
pressions (j40|) and (|4Tj) generalize the formulas (f3Tj) and 
(|32[) . The numerical results for three-copy purification 
are plotted in Fig. [S] This figure confirms the general 
trend that for weak phase fluctuations it is advantageous 
to condition on measurements of p while for strong noise 
it is better to measure the quadrature x [35j . The in- 
termediate strategy where the first balanced homodyne 
detector (BHD) measures x while the second BHD mea- 
sures p does not bring any advantage so we can conclude 
that the optimum strategy that provides maximum re- 
duction of quadrature variance consists of measuring the 
same quadrature (cither x or p) by both homodyne de- 
tectors. 



IV. ASYMPTOTIC LIMIT 

It is possible to find a state to which the purification 
procedure pip converges in the limit of infinite number 
of iterations, k — > 00. This asymptotic state can serve as 
a reference for the performance of the purification, giving 
values that can never be surpassed by finite number of 
iterations. Due to the nature of the result we seek, we 
shall consider the post-selection threshold X in Eq. (fl4|) 
to be zero. That is, given two copies of a noisy state pm,j, 
j = 1,2, the state after a single step of the purification 
is given by the map 

Pout = £|x=0)G°in,l ® Pm,2) 

= Tr 2 [U 12PinA 9 Pi n,2U f 12 \x = 0) 2 (:r = 0|], (42) 

where f 12 denotes a unitary operator describing balanced 
beamsplitter coupling modes 1 and 2 and \x = 0)2(2; = 0| 



8 




FIG. 9: (Color online) Equivalence of the single-step of the 
purification for homodyning (left) and projection onto vac- 
uum (right) as a conditioning measurement. 

is a projector on the eigenstate of the x quadrature. In 
the following we shall exploit the possibility of expressing 
the x eigenstate as 

\x = 0)= lim S r \0), (43) 

r — >oo 

where the operator S r = cxp[r(a 2 — a^ 2 )/2] stands for 
finite squeezing with parameter r and |0) is the vacuum 
state [271 ]. Since the operator SV,i ® S r ,2 commutes with 
the beamsplitter transformation U±2, we can rewrite the 
relation (|4"2"|) as 

p ou t = lim S r £\ ){p' i lr <S> p' i 2r )S}, 

r — ► oc 

= limS r (0| a ^i a ^ n>lr ® / 4 iar ^ a |0) 2 Sj (44) 

with p[ n j r = SJ.pi n jS r and where the purification map 
£|o\ is the same as in Ref. [26fl . This is schematically 
demonstrated in Fig. [9] Note that the inverse squeezing 
operation performed on the output state cancels out with 
the initial squeezing of the next purification step. Thus, 
the iterative purification scheme, using projection onto 
the quadrature eigenvectors in all steps, is equivalent to 
iterative purification with projection onto vacuum states, 
accompanied by squeezing of all initial states and inverse 
squeezing of the purification outcome [l?!], see Fig. [10] 
This allows us to find the asymptotic limit, analogously 
as in Rcf. 26]. Due to this similarity, we shall mention 
here only the most important parts of the derivation and 
the reader can find more details in Ref. [2o| . 

The key realization in finding the limit is that sta- 
tionary points of the purification map £\o){p <8> p), i.e. 
states that remain unchanged by the map, are Gaussian 
states with zero mean value. The initial non-Gaussian 
state will, over the course of purification, converge to 
one of these states. As another important observation, 
let us note that the coefficients oi 1,020 and 002 defined 
as Uij — pij/poo are also unchanged by the purification 
map £|q) and there is a one-to-one correspondence be- 
tween these three coefficients and the set of covariance 
matrices, which in turn completely describe any single- 
mode Gaussian state with zero mean. 

The procedure of finding the limit is the following: 
first, the initial states are described by the Wigner func- 




S I s s s 



\ 1 1 r 

FIG. 10: (Color online) Equivalence of the iterative purifica- 
tion setup for homodyning (left) and projection onto vacuum 
(right) as a conditioning measurements. 

tion 

(45) 

where the vector £ = (x,p) comprises the two conjugate 
quadratures and 

Sr = ( e l ) (46) 

and R(4>) is defined in Eq. (|3"3"|) . The Q-function corre- 
sponding to this Wigner function reads 

Q(x,p) = J $(0)43 ex p (-^ r ^ T ) d< t>> ( 47 ) 

with V r ,<p = (£r,</> + ^/2) _1 , where I stands for the iden- 
tity matrix. Now, with the substitution x = (a+a*)/ v2, 
p = (a — a*)/iy/2 we can use relation ([7]) and find the 
density matrix elements 




where ()^ denotes averaging over the distribution $(^>). 
From the density matrix elements we can obtain the co- 
efficients 011,020 aiml <J 02: G ij — Pij/poo, which are in- 
variant over the course of the purification and, therefore, 
describe the Gaussian limit to which the state converges. 
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It is now a simple matter of finding the inverse relations 



Tn = I - an 
^22 = 1 — en 

T21 = Ti2 = 



O"02 + CT20 
<702 + C20 



V2 

O"02 — 020 

iV2 



(49) 



After some algebra we find that the formulas (|48|) and 
(|4*9")l can be simplified to 



(ry,0 V|r r ,0|)0 
(x/ir^D^ 



(50) 



The final covariance matrix is given by E = S'~ 1 (r _1 — 
I/2)S~ 1 . Due to the symmetric nature of the phase 
noise, the matrix E is diagonal and fully specified by val- 
ues of variances of the squeezed and antisqueezed quadra- 
tures, V^ii m and VpJim- We can directly find the output 
variances by taking the limit r — > oo, obtaining 



V, 



V, 



p,lim 



v x v p 



(AW 



where 



A = (V x cos 2 



V p sin 



-1/2 



(51) 



(52) 



Note that Kc.iim^pjiim = V X V P so the purification proce- 
dure asymptotically suppresses all added noise and the 
asymptotic Gaussian state has the same purity V = 
l/(2y^V x Vp) as the initial state before de-phasing. This, 
however, holds only for idealized conditioning on x = 0. 
For finite acceptance window, X > 0, the purity of the 
asymptotic distilled state would generally be lower than 
the purity of the initial state. 

It is possible to modify the method to incorporate 
other scenarios discussed in previous sections. Find- 
ing the limit state for iterations incorporating condition- 
ing measurements along an arbitrary quadrature q{9) is 
straightforward, requiring only the use of an alternative 
squeezing operation described by S r ,0 = R(0)S r R T (8). 
Adapting the calculations to incorporate for imperfect 
detectors is a bit more tricky. Inefficient detection with 
efficiency rj can be simulated by a beamsplitter with 
transmittance rj followed by a perfect detector. The at- 
tenuation can be represented by a map, M(p), acting on 
the density matrix and transforming the corresponding 
covariance matrix of the Gaussian state as, 



1 — 77 

E - 77E + -^ij. 

The inverse operation, A4~ 1 (p) acting as 

E^i(E-1^7) 
77 2 



(53) 



(54) 





FIG. 11: (Color online) Performance of the iterative purifica- 
tion protocol (green line), compared to the initial phase dif- 
fused state (blue dashed line), asymptotic limit (red circles) 
and state without phase diffusion (black plus signs). The 
arrows indicate increasing number of iterations. The values 
shown are as if verified by an imperfect detector. The param- 
eters were V x = 0.2, V p = 2, 77 = 0.85 and X = 0. 



is, of course, nonphysical. We can, however, formally use 
it and in analogy with Eq. (|4~4"|) represent a single step of 
the purification procedure by the relation 

Pout = M- 1 ( lim S r sSl) , (55) 

where the e is given by 

S = £| ) (5+ ^(pin^SrA ® S^O'm.a)^) -(56) 

Now again, starting from the modified input states 
SlAi(p)S r we can find the asymptotic limit for the map 
£|o> and with the use of relation (|55)) it can be trans- 
formed into the desired result. However, if we are inter- 
ested in results that can be observed by imperfect detec- 
tion, the final operation Mr 1 is not necessary. 

Some numerical results are shown in Fig. 1111 It can 
be seen that the iterations converge to the limit given 
by Eq. (|51|) in both the resulting variance and purity. 
Moreover, the purity of the limit is independent of the 
actual amount of initial phase fluctuations and equals to 
the purity of the original state (modified by imperfect de- 
tectors), although it is less when the states are compared 
directly, without the imperfect verification measurement. 

Figures [T2l and [T3l illustrate the behavior of the asymp- 
totic limit relative to the angle 9 of the conditioning mea- 
surement. Although all choices of 9 lead to a state puri- 
fied to some extent, the two most interesting choices are 
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FIG. 12: (Color online) Variance of the asymptotic limit state 
relative to amount of phase fluctuations, a and angle of con- 
ditioning detection, 6. The values shown are as if verified by 
an imperfect detector. The parameters were V x = 0.2, V p = 2 
and rj — 0.85. 




FIG. 13: (Color online) Purity of the asymptotic limit state 
relative to amount of phase fluctuations, a and angle of con- 
ditioning detection, 9. The values shown are as if verified by 
an imperfect detector. The parameters were V x = 0.2, V p = 2 
and rj = 0.85. 

9 = 0, corresponding to the measurement of the quadra- 
ture operator x and 9 = n/2, corresponding to the mea- 
surement of p. Only these two choices can lead to a purity 



of the final state to be maximal, that is, equal to the pu- 
rity of the initial state (modified by detection efficiency). 
This is not very surprising, since these two choices coin- 
cide with the basis in which the covariance matrix of the 
initial state is diagonal. We can see that for weak phase 
noise (small a) it is advantageous to employ condition- 
ing on measurements of the p quadrature while for strong 
phase fluctuations the x based post-selection is optimal. 
V. SUMMARY 

We have provided a detailed analysis of multiple 
copy purification protocols for phase diffused single-mode 
squeezed states. We performed most calculations in the 
Fock basis, which allowed for efficient numerical treat- 
ment. We have also discussed the performance of col- 
lective purification, employing more than two copies of 
the initial state in one step. The collective purifica- 
tion/distillation uses an arbitrary number of copies to 
improve the purity and non-classicality of the state, and 
therefore was found to be a generalization of iterative 
purification and distillation which requires 2 fe copies. 
In particular we have shown that collective distillation 
based on three copies of phase diffused squeezed states 
is possible, and already provides a significantly improved 
degree of squeezing when compared with two-copy dis- 
tillation. We have also investigated the convergence 
of the iterative purification protocol and derived semi- 
analytical expressions for variances of the asymptotic 
Gaussian state. Finally, let us note that in the spirit 
of rcf. [lH, the results presented here can be straight- 
forwardly extended to describe the distillation of phase 
diffused two-mode squeezed states. 
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